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Abstract: A supersymmetric extension of the nonlinear 0(3) sigma model in two spacetime 
dimensions is investigated by means of Monte Cai^lo simulations. We argue that it is impossible 
to construct a lattice action that implements both the 0(3) symmetry as well as at least one su- 
persymmetry exactly at finite lattice spacing. It is shown by explicit calculations that previously 
proposed discretizations fail to reproduce the exact symmetries of the target manifold in the con- 
tinuum limit. We provide an alternative lattice action with exact 0(3) symmetry and compare 
two approaches based on different derivative operators. Using the nonlocal SLAC derivative for 
the quenched model on moderately sized lattices we extract the value cj(2,iio) = 1.2604(13) for 
the step scaling function at uq = 1.0595, to be compared with the exact value 1.261210. For 
the supersymmetric model with SLAC derivative the discrete chiral symmetry is maintained but 
we encounter strong sign fluctuations, rendering large lattice simulations ineffective. By applying 
the Wilson prescription, supersymmetry and chiral symmetry are broken explicitly at finite lattice 
spacing, though there is clear evidence that both are restored in the continuum limit by fine tuning 
of a single mass parameter. 
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1 Introduction 



Despite its great success in the description of high energy experiments, the Standard Model of 
particle physics still faces serious problems, especially with regard to cosmological observations, 
where yet undiscovered sources of gravitation are proposed to correct the seemingly unbalanced 
equation of observed matter on the one hand and gravitational attraction on the other hand [1]. 
Another interesting question concerns the existence of far more matter than antimatter in our uni- 
verse, and this fact is connected to the strong CP problem (see [2] and references listed therein). 
Coined as the hierarchy problem is the question of scales that are involved in the fundamental the- 
ory of particle physics. To address these open issues, extensions of the Standard Model have been 
proposed, among them supersymmetrically extended versions, in particular the Minimal Supersym- 
metric Standard Model (MSSM). It is of great interest to examine the nonperturbative properties 
of this theory, but the utilization of the only ab initio method for this purpose, namely a spacetime 
lattice simulation, is notoriously difficult because of the many different scales that are involved. 
In the past, two-dimensional nonlinear sigma models have been applied successfully to model non- 
perturbative properties of four-dimensional strongly coupled pure gauge theories [3]. Similarly 
one may employ the supersymmetrized version of the nonlinear sigma model to effectively de- 
scribe super- Yang-Mills theories with a strongly interacting fermionic sector. A first construction 
of the supersymmetric nonlinear sigma model with O(A^) target manifold is due to E. Witten [4] 
and P. Di Vecchia and S. Fenaia [5], followed by a number of papers that established analytical 
properties, including asymptotic safety, spontaneous breaking of chiral symmetry and the dynami- 
cal generation of particle masses [6-8]. Further studies are particularly concerned with the special 
case of the 0(3) model, which admits an extended J\f=2 supersymmetry algebra since its target 
manifold is Kahler [9]. 

This enhanced supersymmetry has recently raised some interest with regard to the endeavor to study 
supersymmetric theories by numerical simulations. Because it is an extension of the Poincare sym- 
metry, any lattice discretization of spacetime also breaks supersymmetry. This can be traced back 
to the failure of the Leibniz rule on the lattice [10]. It is therefore necessary to take great care in 
restoring the full symmetry in the continuum limit. In general, one needs to introduce appropriate 
counter terms for each relevant supersymmetry-breaking operator that must each be fine tuned such 
that the desired continuum theory is reached [11]. Depending on the number of parameters, this be- 
comes practically impossible very fast due to hmited computer time. In addition, much information 
about the theory is needed prior to numerical investigations. For that reason alternative approaches 
have been proposed in order to reduce the number of fine tuning parameters or even render the 
fine tuning procedure obsolete [12]. One such approach for theories with extended supersymmetry 
aims at the construction of a nilpotent charge composed from the supercharges such that both the 
continuum as well as the discretized model ai^e invariant under this charge. It is then expected that, 
by preserving a part of the symmetry, invariance under the full supersymmetry is restored automat- 
ically in the continuum limit, without the need for fine tuning. This procedure was applied to the 
supersymmetric 0(3) model [13, 14] and the authors conclude that the supersymmetric Ward iden- 
tities are indeed fulfilled in the continuum limit, hinting at a restoration of the full supersymmetry. 
However, the lattice discretization constructed in this way breaks the 0(3) symmetry of the target 
space explicitly at finite lattice spacing and no attempt was made to show that it is restored in the 
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continuum limit. It is therefore unclear whether this discretization actually describes the correct 
model in the continuum and we will address this question in much detail. 

Expensive dynamical fermion simulations are needed in order to treat the superpartners on equal 
footing and we wish to utilize the SLAC prescription for the fermionic derivative, which performed 
very well in previously examined Wess-Zumino models [15], showing greatly reduced lattice ai^- 
tifacts in comparison to e.g. the Wilson derivative. Furthermore, the SLAC derivative allows for 
a discretized action that does not break the discrete chiral symmetry explicitly. But this comes at 
the cost of having a nonlocal derivative operator and we will explicitly check whether the SLAC 
derivative is applicable in the present model. In order to compare with an alternative, the model 
will also be studied by simulations based on the Wilson derivative, which is less affected by the 
sign problem and enables us to investigate larger lattices. 

The article is organized as follows: We introduce the supersymmetric version of the nonlinear 0(3) 
model along with a derivation of its stereographic projection in the superfield formalism. We then 
discuss the extended supersymmetry of the model and express the related transformations in terms 
of constrained fields. Based on the continuum theory two alternative discretizations in Euclidean 
space are presented which resolve the field constraints - one using stereographically projected 
fields and the other using elements of the orthogonal group S0(3). Thus providing the basic setup 
of our lattice discretization, we go on to argue that the SLAC derivative is indeed applicable here by 
performing a high precision analysis of the step scaling function in the quenched model. Towards 
the end of the first part we finally answer the question whether a lattice theory exists that admits 
both the 0(3) symmetry and at least one supersymmetry. 

In the second part, which is primarily dedicated to extensive numerical investigations, we discuss 
the missing target space symmetry in a previously given lattice discretization by Catterall et al. 
[13] and proceed with results from our own calculations. We examine the discrete chiral symmetry 
as well as supersymmetry by presenting results for the chiral condensate and the masses of the ele- 
mentary excitations computed from ensembles obtained with both SLAC and Wilson derivative. In 
the second case special emphasis is put on the fine tuning of divergent operators in order to anive 
at a supersymmetric continuum limit. Before coming to our conclusions, we present our results for 
a simple Ward identity based on the expectation value of the bosonic action. 

2 Symmetries and possible discretizations of the O(A^) nonlinear sigma model 
2.1 Constrained and unconstrained formulations 

The supersymmetric extension of the nonlinear O(A^) sigma model in two-dimensional Euclidean 
spacetime can be formulated in terms of a real superfield ^ , 

^ = n + i6'tl) + ^e9f, (2.1) 

subject to the constraint = 1, where n and / ai^e A^-tuples of real scalar- fields and i/j is an 
A^-tuple of Majorana fields. We shall refer to the elements of a tupel as "flavors". The constraint in 
superspace is equivalent to the following constraints for the component fields, 

= 1, nip = and nf = ^ipil^ . (2.2) 
'Conventions concerning the Gamma matrices and the Fierz relations are explained in appendix A. 
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The 0(A^)-invariant Lagrangian density is defined in terms of the covariant derivatives = 

C = ^-mD^\ss = ^ {d^nd'^n + iV'^V' - /') ■ (2-3) 

The equation of motion for the auxiliary field implies that / and n are parallel, / = ^{'il)ip)n, and 
the resulting on-shell Lagrangian density contains a 4-Fermi term, 

/: = ^ [^^n^^'n + + i (V^^)^) . (2.4) 

The action and the constraints ai^e both invariant under global 0{N) "flavor" transformations and 
by construction they ai^e also invariant under the = 1 supersymmetry transformations 

5n = iei/^, Sip = [0 + ^ipip) ne . (2.5) 

Besides flavor symmetry and supersymmetry the classical theory admits a further Z2-symmetry 
generated by the chiral transformation -0 — However, quantum fluctuations dynamically 
generate a mass term and hence induce spontaneous breaking of the chiral Z2-symmetry. 
Explicit constraints for the fields are sometimes difficult to handle, e.g. in numerical simulations. 
Therefore, it is useful to construct a formulation of the model in terms of an unconstrained real 
superfield U{x, 9) = u{x) + i^A(x) + ^99g{x), which is an {N — l)-tupel. It is related to the 
superfield $ by a stereographic projection in superspace 



$1 \ 1 



(2.6) 



I *_Ly i + i72 1 2(7 / ■ 

The decomposition of the projection into bosonic and fermionic fields reads: 

n± = 2pu, tp± = 2pX — ip^ (uX)u with p= ;r . (2.7) 

1 + 

The expressions for the remaining components ni and ipi can be determined either from (2.6) or 

from (2.7) and the constraints -n? = I and mp = 0. The inverse transformation in superspace 
reads U = (1 + $1) and leads to 

u = ^nj_, X = -^ip±- ^ipinj_ with p= ^ . (2.8) 

Applying the stereographic projection (2.7), the on-shell Lagrangian density can be written in terms 
of the unconstrained fields as 

C = ^p^ {d^ud^'u + iA^A + 4ip {Xu)-f^'{d^uX) + p2(AA)2) . (2.9) 

The action is invariant under the supersymmetry transformations 

6u = ieX, 6X= {$ + ipXX)ue-2ip{Xu)Xe. (2.10) 
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2.2 Extended supersymmetry of the 0(3)-model 

The fields u and A are unconstrained, but their target manifold has a non-trivial metric. In fact, 
for = 3 the target manifold is Kahler and the corresponding potential can be written in terms of 
the complex field u = ui + iu2 as K{u, u) = ln(l + uu). As pointed out by Zumino a nonlinear 
sigma model whose tai^get manifold is Kahler, possesses an AA=2-supersymmetric extension [9]. 
Hence the 0(3) admits an additional supersymmetry besides (2.10) and (2.5). 
In order to derive an explicit expression of this symmetry, we investigate a general ansatz in terms 
of the unconstrained fields (6u = ie{Ai)X, etc.), and determine the constraints for the matrices 
Aj, etc., which follow from the supersymmetry algebra and the invariance of the action^. This 
approach yields the second pair of symmetry transformations as 

6u = a2eX (2.11) 
6\ = i(j2 {(^u — ip {\\)u + 2ip (Ait)A) e . 

Both supersymmetries (2.10,2.11) can also be obtained by deriving the complex supersymmetry 
from the Kahler potential, cf. [9], and decomposing the complex fields and complex transforma- 
tion parameters into real ones. 

Applying (2.7) and (2.8), we can also express these transformations in terms of the constrained 
fields in order to determine the second supersymmetry of (2.4). One finds the simple transforma- 
tions 

5n= in X ei/; (2.12) 
6ijj = —n X d^n^^e — leip x i/) , 

where a x b denotes the vector product of a and b. A proof that the action (2.4) is invariant 
under these transformations is given in appendix B. The two on-shell supersymmetries (2.5,2.12) 
are generated by the supercharges 

Qi = l^fd^ml^ , Qu = -ij 7''7°(« x d^n)xl^ . (2.13) 
These results are in agreement with the supercurrents constructed in [4]. 

2.3 Discretization and constraints 

So far, sigma models in the continuum have been considered. In order to investigate the corre- 
sponding lattice models one should try to discretize it in a way that maintains as many symmetries 
of the continuum theory as possible. This is difficult with respect to supersymmetry, as will be 
discussed in section 2.5. But also the flavor symmetry must be treated with care: If one starts with 
an unconstrained formulation of the model and tries to discretize it in a straightforward way, one 
generically breaks the 0{N) symmetry. For the 0(3) model we will illustrate this in more detail 
in section 3.1. Simulations demonstrate clearly that even in the continuum limit the symmetry is 

^An example of this approach can be found in [16], where it is applied to the Wess-Zumino model. 



-5- 



not restored in such cases. In order to avoid this problem, we start with a formulation in terms of 
constrained fields, whose discretization is manifestly 0(A^)-invariant: 

5[n,^] = ^ {nlK^yny + i^^M^yP^^ + liil^Ay^y)^). (2.14) 

The subscripts x, y denote lattice sites, while a, /3 are spinor indices. The lattice derivatives Kxy 
and Mxy are proportional to the identity in flavor space. The constraints UxUx = 1 and n^il^x = 
must be fulfilled at each lattice point x and they ai^e implemented as delta-functions in the path 
integral measure. This causes some difficulties in numerical simulations. One can cope with this 
problem by either applying the stereographic projection or by introducing group valued dynamical 
variables. 

2.3.1 Stereographic projection and measure of path integration 

The stereographic projection (2.7) resolves both constraints and leads to an unconstrained but yet 
0(A^)-symmetric lattice formulation: 

S[u, A] = + S2F + 5'4F, with 

Y^Yl ^PxUlKxyUyPy + Px{l - ul)Kxy{l - Ul)py, 

x,y 
x,y;a,l3 

Saf = -^^.pti^xK)'^ ■ 

^ X 

(2.15) 

The change from the constrained fields (n, to the unconstrained fields (u, A) yields a non-trivial 
Jacobian which is computed in appendix C. The result is^: 

n dn, ^(^2 _ X)b{rxi\^x)b(ri^l) J] ^A^ (l + tx^)^-^ . (2.16) 

The four-fermion interaction can be eliminated by employing the usual Hubbard-Stratonovich 
transformation [17], which introduces an auxiliary bosonic field a: 

S\u, A] = Sb + 52F + ^ ^ {ol + ^xuxpI ^xK) (2. 17) 

2.3.2 Coset formulation 

The constrained field n propagates on the unit sphere S^' with A^' = A^ — 1 which can be viewed as 
coset space SO(A^)/SO(A^'). To relate the constrained field to an element of the orthogonal group 
SO(A^) we supplement at each lattice site the unit vector n by orthonormal vectors ei, . . . , e^v' 

^One obtains the same result if one does not eliminate the auxiliary field / at the beginning but projects it in accor- 
dance to (2.6) and integrates out the unconstrained auxiliary field g afterwards. Following this approach, the superdeter- 
minant yields only a trival factor, while the density 1/ p^^^ stems from integrating w.r.t. g. 



2p^uu^)^MSy^{p 



2p'uu^) + 4 {p'u)^M^^ {p'u^)^ 
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such that {n, ei, . . . , e^r/} forms an oriented orthonormal basis of R^. The A^-component Majo- 
rana spinor orthogonal to n is a linear combination of the Cj, i.e. 

^ = XieiH ^Xn'^n'- (2.18) 

These N basis vectors may be viewed as columns of a rotation matrix R = (n, ei, . . . , e^r/) in 
N dimensions. Therefore the x-dependent orthonormal basis is given by a spacetime-dependent 
rotation R G SO(A^) acting on a constant (x-independent) orthonormal frame {no, gi, ■ ■ ■ ,gN'}' 

n{x) = R{x)no, ei{x) = R{x)gi , (2.19) 

and the path integral in terms of new variables R, x ^rid a reads 

Z = y^Pi^Po-Pxe-^B-Sp (2.20) 

with actions 

-Sb = ^ ^ {nlRlK^yRyno + a^S^yay) , 5"? = ^ ^ xlQxyXy ■ (2-21) 

The action for the fermions contains the real and antisymmetric matrix 

Qxy,ij = gj R^CMxyRygj + SxySijCcTx ■ (2.22) 

The base vectors {g^} (or {ej}) are not unique since any rotation in the plane orthogonal to no 
transforms an admissible set of base vectors into another admissible set. More precisely, 

R — y R' = RS, Xi — >Xi = XkSki, with 5no = no, Sij = gJSgj 

are local SO(A^') symmetries of the action (2.21). The measure of the path integral is not affected 
by the change of the dynamical fields from n to R, i.e. a distribution of R according to the Haar 
measure on SO(A^) will lead to a uniform distribution of n = Ruq on the sphere. This renders the 
structure of the lattice action rather simple. The price we pay is a local SO(A^') symmetry in the 
choice of the basis vectors e^. 

2.3.3 Fermion determinant 

The treatment of the fermion operator is very similar for stereographically projected or group val- 
ued variables and we will only depict the second case here. Performing the Grassmannian integral 
in (2.20) leads to the bosonic path integral 

Z = ypflPcj signPf(g)e~^B+in|Pf(Q)| ^ (2.23) 

and simulations are performed in the sign quenched ensemble. This means that the configurations 
are sampled according to the probability distribution defined by the exponential function. The sign 
of the Pfaffian is handled by reweighting. In this formulation the invariance of the path integral 
under the local SO(A^') transformation is obvious: The bosonic action and Haar measure are both 
invariant under the change of variables R — RS with S'no = no. The operator Q^y transforms 
as SjQxySy, which leaves the Pfaffian invariant since S = ^xSx cancels in the general relation 
Pi{S'^QS) = det(5)Pf(Q). In the simulations the effective fermionic action is rewritten ac- 
cording to In |Pf(Q)| = i Indet(Q) and the hybrid Monte Carlo algorithm is used, such that the 
Hamiltonian evolution of the group valued R field is similar to the quenched case. 
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2.4 Lattice derivatives 

Now we specify the lattice derivative for the fermions Mxy and bosons K^y in the 0{N) invariant 
formulation. We will use two different implementations, the SLAC derivative and the Wilson 
derivative, so that we can compare results. 

2.4.1 SLAC derivative 

With the help of discrete Fourier transformation on a finite A't x N^; lattice A with N lattice sites 
one may motivate the SLAC derivative [18], 

-j= E ^p,-'''np) = E f(y) E ^p,-''^''-'^) = E ^ (2-24) 

^ peA* yeA pgA* j/GA 

where p is from the dual Lattice A*. This non-local derivative has been proven to be useful in the 
context of Wess-Zumino models [15] and supersymmetric quantum mechanics [19, 20]. Because 
the formulations based on the SLAC derivative show only small lattice artifacts and do not break 
the ^2 chiral symmetry explicitly, we use this derivative for the supersymmetric nonlinear sigma 
model: 

K.y = -Y. i^f^nly, Mxy = (T-^r^),^ . (2.25) 
The SLAC -derivative leads to lattice models without doublers. 

2.4.2 Wilson derivative 

Another suitable choice that avoids fermion doublers is the Wilson derivative which introduces an 
additional momentum-dependent mass term that vanishes in the naive continuum limit, 

M^y^ = 7f{d^r).y + ^""^Y^^y ' (2-26) 

where 5^^™ is the symmetric lattice derivative, a the lattice spacing and A^y the lattice Laplacian. 
We choose the bosonic derivative in a particular form that has demonstrated superior results in 
Wess-Zumino models [15], 

^^. = -E(^r)', + (y^^.)'- (2-2V) 

A* 

The Wilson prescription breaks chiral symmetry explicitly, but is ultralocal in contrast to the SLAC 
derivative. 

2.4.3 Universality and the SLAC derivative 

It is not obvious that the SLAC derivative can be used for models with curved target space. To 
justify its use we will first study the purely bosonic 0(3)-model with this derivative and test if 
the known continuum result with scaling of the finite volume mass gap m can be reproduced. A 
quantity accessible at finite volumes is the step scaling function introduced by Liischer, Weisz and 
Wolff [21]. At finite spatial volume L = aNf^, the renormalization group invariant variable m{Nfi)L 



-8- 



can be computed at large temporal extent. On every fi- 
nite lattice the step scaling function S is then determined 
according to 

S(2, u, N-') = m(2iV,)2iL(jvoL=. > (2-28) 

where the condition m{N^)L = u determines the bare 
coupling g^"^ that is used on both spatial volumes A'^s and 
2iVs. 

One has to measure the finite volume mass gap m twice 
in order to determine S at a given N^. First one has to 
tune the bare coupling until mL is u on the small lattice 
with spatial sites. Then one goes to the larger lattice 
with 2A'^5 sites and measures m{2N^)2L^ with the previ- 
ously determined bare coupling. The lattice step scaling 
function is expected to have a universal continuum limit 

(j{2,u) = S(2,u,0) = lim S(2,n,iV-^) 

and the high-precision results in [22, 23] demonstrate this universal behaviour. Here, the step 
scaling function is computed using the SLAC derivative. Following [21] we use uq = 1.0595 
which, according to [24], leads to the continuum value 

C7(2, no) = 1.261210. (2.29) 

Using the SLAC derivative the regularized action is given by 

S = -^Y1 nlRl{df^'^^)%Ryno, (2.30) 

x,y 

where R is the group valued dynamical field. Since the action is not given by nearest neighbour 
interactions, a cluster algorithm is not applicable and the hybrid Monte-Carlo algorithm will be 
used instead. The momenta of the 'Hamiltonian' used for the hybrid Monte-Carlo algorithm aie 
elements of the Lie algebra so(3). 
The 0(3) invariant correlator is naturally defined as 

C(t) = iV-2^(n(i,,)n(o,,)> (2.31) 

x,y 

and the mass is extracted via a fit to 

C(t) oc cosh (ma (t-iVt/2)) (2.32) 

on a logarithmic scale, so that contributions from t's near are taken into account, where the 
contributions of the higher excited states are suppressed. However, reliable high-precision results 
can only be obtained if systematic errors are under control. To actually see the contribution of 
higher excited states, one considers the extracted mass as a function of to where the fit is performed 
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1.26587(35) 
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20 


1.26416(70) 


1.52870(12) 



Table 1. Value of the step scaling func- 
tion Yj{2,uq,N^^) after interpolating to 
uq = 1.0595 for different spatial lattice 
sizes and corresponding g^^ at which 
m{N^)L ~ uq is reached. 
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over a range t G [to, — to] for fixed Ni and A'^g. This effect is analysed for A'^s G {6, 12} and 
A't = 6iVs for a coupling g~'^ = 1.309, which is quite close to the point of interest (2.29), with 
extremely large statistics of about 5 • 10^ configurations, distributed over 1 000 replica. The results 
depicted in Fig. 1 show that for larger lattices the contribution of higher excited states is well below 
the usual statistical accuracy that is used for most of the results given below. For small lattices there 
ai^e two competing effects: For small to the contribution of higher states is well visible, whereas for 
large to the well-known fluctuations arising from the non-locality of the SLAC derivative begin to 
grow. Therefore the optimal choice leading to systematic eiTors of the same order as the statistical 
ones, is given by to = -^s and will be used in the following. The second possible systematic en^or 
is due to the finite A't. For A't too small a thermal contribution to the mass gap will show up. This 
effect is investigated by keeping to = -^s fixed but varying A't, see Fig. 2 for A'^^ G {6, 12}. For 
the smaller lattice the contributions at small A't are more pronounced and become negligible for 
A't > 6Nfi while for the larger lattice A't > 5A'^s is sufficient. To suppress these systematic eiTors 
A't = SA'^s (A'^t = 6Nfi) is used on the smaller (larger) lattice of each step scaling computation. 



1.0630 



1.0625 



mL 



1.0620 



1.0615 



1.0610 




1.286 



mL 




Figure 1. Mass gap extracted from a logarithmic cosh fit of the correlator (2.31) in a range t G [tQ,Nt — to] 
for A^s = 6 (left panel) and = 12 (right panel) at coupling g^^ = 1.309 for fixed A't = GN^- The shaded 
area denotes the usual accuracy of results at other g^^. 
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Figure 2. Mass gap extracted from a logarithmic cosh fit of the correlator (2.3 1) in a range t G [A^s, A't — A^s 
for Ns^6 (left panel) and N, 12 (right panel) at coupling ^ = 1.309 for different Nt/N,. 



-10- 



With the systematic errors under control it is possible to extrapolate the lattice step scaling function 
to the continuum limit. The bare coupling is tuned such that u = m{Ns)L is near uq = 1.0595 on 
lattices with sizes G {4, 6, 8, 10, 12, 16, 20}. The corresponding S(2, u, N~'^) are plotted over 
u for a subset of these A'^s in Fig. 3 (left panel), and a linear interpolation based on seven different 
coupling g~'^ allows for the extraction of S(2, uq, N^^) at the point uq = 1.0595. The explicit 
values are collected in Tab. 1 . With Symanzik's theory of lattice artefacts it has been argued in 
[22] that finite a corrections are of order 0(a^(lna)^) and appear nearly linear for a large range 
of computationally accessible lattice sizes [25]. For that reason an extrapolation to a = based on 
the formula ^ ^ 

j:{2,uo,N~') = a{2,uo)+A(^-^^ (^In^^ (2.33) 

is used. The results for A't = 4 have been omitted because of the large systematic errors introduced 
by the fluctuations arising from the SLAC derivative for large lattice spacings. The extrapolation 
is shown in Fig. 3 (right panel) and a value of cr(2, uq) = 1.2604(13) is extracted. This value is in 
complete agreement with the exactly known result in the continuum limit, see eq. (2.29). Therefore 
a discretisation of the (bosonic) 0(3) nonlinear sigma model with the SLAC derivative is feasible 
and may also be used for the supersymmetric model. 
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Figure 3. Left panel: Step scaling function for lattices with spatial size A^s G {4, 6, 8, 12}. Shaded regions 
indicate error bounds of a linear interpolation. Right panel: Continuum limit of the step scaling function for 
Uq = 1.0595. The shaded area indicates the error bounds of a fit according to Eq. (2.33), where the value 
for = 0.25 has been omitted. The black dot marks the continuum value given in Eq. (2.29). 



2.5 Supersymmetries on the lattice 

Is it possible to discretize the supersymmetric nonlinear 0(3) -model such that its characteristic 
symmetries are maintained? We have seen that it is straightforward to find a manifestly 0(3)- 
symmetric discretization. The treatment of supersymmetry on the lattice is more difficult, because 
it is an extension of the Poincare symmetry, which is broken by any discretization of spacetime. 
However, there are some approaches to maintain at least a part of the supersymmetry on the lattice. 
One of these approaches relies on a twisting of the supersymmetry in a way that provides a nilpotent 
supercharge Q which is used to find a Q-exact formulation S = QA of the lattice action. By 
construction this action is invariant under supersymmetries generated by Q. This approach has 
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been applied to the 0(3)-NLSM in [13, 14]. However, the authors employ a formulation of the 
model in terms of unconstrained fields and its discretization breaks the 0(3) symmetry in such a 
way that it is not restored in the continuum limit. We will demonstrate this explicitly in section 3.1. 
We must thus conclude that the continuum limits of these lattice constructions cannot be identified 
with the two-dimensional non-lineai^ 0(3) sigma model. 

Are there other ways to find a partly supersymmetric but still 0(3) symmetric discretization? A 
symmetry of the model has to be a symmetry of the action (2.3), but is also has to be compatible 
with the constraints n? = 1 and mp = 0. Any supersymmetry has to be a combination of the 
transformations given in (2.5) and (2.12). If we consider the discretization of these, we notice that 
the first transformation (2.5) breaks the constraint nip = on the lattice, because 

5i(n,C) = + E ^-D^y'^y^^ + K^.-V'x>^e° ^^=^ ^ n^D^^UyS^ (2.34) 

yeA j/GA 

does not vanish for arbitrary n^., no matter which derivative Dxy we use^. In contrast, the second 
transformation respects the constraints at each point: 

Siiirixtp^) = [{n^ X £-0^)^" - ^n^in^ x D^yUye'^) - m^ieip^ x i/?^) = 

s/gA 

5ii{nl) = 2in^(n^. x e^^) = . (2.35) 

We conclude that no nontrivial combination of the two transformations 5i and 5ii can be a symmetry 
of the lattice theory, since the second transformation cannot restore the violation of the constraints 
caused by the first one. The second transformation on its own, however, can also not be a symme- 
try of the action because of {Q^ , Q}^} = 2^7^^ 9^. The superalgebra furthermore tells us that an 
approach based on a nilpotent supercharge is not possible, either, because such a nilpotent charge 
would have to be a combination of both supercharges and would hence violate the constraints. 
Could we circumvent this restriction by improving the discretiziation in some way? Comparing our 
formulation with the one investigated in [13, 14], we see that the latter one contains an additional 
topological term. However, such a term does not affect the supersymmetry transformations (2.5) 
and (2.12) and hence cannot solve the problem. From a systematic point of view, there are only 
two modifications possible which are compatible with an 0(3)-invariant continuum limit. The first 
possibility is to modify the terms that are already present in the action. For example, one could 
introduce non-local interaction terms like Cxyzw{'4>xtpy)iipzipw) instead of l]^(i/'z'0s)^ 

[26]. The second possibility could be an inclusion of additional terms in the lattice action which 
vanish in the continuum limit. Any change of the action, however, does not have an impact on the 
constraints and hence cannot prevent their breaking. A modification of the constraints, by contrast, 
would directly alter the geometry of the target manifold and is thus no alternative. It follows that 
an improvement of the discretization could only maintain a part of supersymmetry by rendering 
the lattice action invariant under the second transformations. But this is not possible due to the 
structure of the superalgebra. 

Even though these arguments were developed for a specific choice of coordinates, they also hold 
true for any reparametrization {n,tp) — > {n',tp'), because such a transformation is a bijective 

''Actually, it is also not a symmetry of the discretized action, but the breaking of the constraints is more severe. 
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mapping between field values at a certain point x in spacetime, which commutes with discretiza- 
tion. As a consequence, one would observe the same pattern of symmetry breaking as depicted in 
(2.34) and (2.35). The single ambiguity that could arise here from the discretized derivative of the 
bosonic field is irrelevant since the presented arguments do not depend on the details of the lattice 
derivatives. 

We conclude that it is just not possible to construct a discretization of the nonlinear 0(3) model 
which maintains 0(3) invariance as well as an exact supersymmetry. From this point of view, the 
symmetry breaking in the ansatz [13, 14] was inevitable. In this article we will work with the lattice 
formulation introduced in the previous section. It maintains the important 0(3) symmetry of the 
model while it breaks all supersymmetries. The latter should be restored in the continuum limit. 

3 Numerical results 

Before we present the results which are obtained by simulations for the 0(3)-symmetric lattice 
models described above, we want to demonstrate how a violation of the 0(3) symmetry arises in 
generic lattice formulations given in terms of unconstrained fields. We shall see that the symmetry 
of the continuum model is not recovered in the continuum limit. 

3.1 Drawbacks of a Q-exact lattice formulation 

This problem can already be illustrated by investigating the bosonic 0(3)-NLSM. If we use a lattice 
derivative that is based on nearest neighbour interactions, the discretization of the bosonic part of 
(2.9) reads 

S = ply\Ux - Uy\'^ ■ (3.1) 

{xy) 

At this point, an ambiguity in defining the density p^y arises. It must interpolate between and 
p^, which coincide in the continuum limit. Naively, many interpolations are feasible, e.g. one 
could use the arithmetic mean, p1.y = i(p^ + p^), or the geometric mean = pxPy In order 
to determine the appropriate interpolation one can apply the stereographic projection on the mani- 
festly 0(3) invariant discretization in terms of the constrained n-field and finds that p1.y has to be 
implemented as the geometric mean. 

To analyze a possible symmetry breaking in the formulation based on the arithmetic mean, sim- 
ulations for both lattice prescriptions have been performed^ and it was tested if (n) = for 
n = A^~^ Ylix ''^x accordance with the global 0(3) symmetry, which cannot be spontaneously 
broken in two spacetime dimensions [27]. 

Fig. 4 shows that the discretization based on the geometric mean indeed yields 0(3) symmetric 
results whereas for the arithmetic mean only a 0(2) symmetry around the 1-axis remains. This 
is a direct consequence of the global U{1) symmetry u — )• e^'^ti, where u = ui + in2. In the 
naive continuum limit both prescriptions are expected to coincide. In order to investigate this issue, 
simulations of the model based on the arithmetic mean have been earned out with different lattice 
sizes N and at different couplings g~'^. A restoration of the 0(3) symmetry implies a vanishing 

""In the hybrid MC algorithm the factor of the path integral measure for the bosonic model, see (C.6), is absorbed 
in the action. 
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Figure 4. Scatter plot (projected to the ni-n^ plane) of the averaged field n for a lattice discretisation based 
on the arithmetic mean (left panel) and geometric mean (right panel) at g^^ = 1 and lattice size = 10 x 10. 



(na) in the continuum limit. For a wide range of couplings (^3) is independent of the lattice sizes 
such that the vacuum expectation value is assumed to be free of finite size effects (see Fig. 5, left 
panel). By fitting the correlator 

Cjiit) = N~^ReY,{uit,y)U^o,x)) (3.2) 
x,y 

to cosh(ma(t — Ni/2)) the mass mL = maN^ measured in units of the physical box length can be 
extracted^. The analysis of (^3) at fixed physical box size mL in the continuum limit (see Fig. 5, 
right panel) clearly shows that {n^) grows to a value close to 1 for small lattice spacings, i.e. for 
large N and large g~^. Therefore it is impossible to reach the coiTcct 0(3) symmetric continuum 
limit for a regularization based on the arithmetic mean. 




Figure 5. The value of (713) as indicator for a broken global 0(3) symmetry for three lattice volumes plotted 
over the couphng g^^ (left panel) and physical box size (right panel). 

*Due to the U(l) symmetry {u) = 0, and Ci.{t) is the connected 2-point function. 
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Having encountered these problems in the simple bosonic case, we ought to be very careful with 
a discretization of the superymmetric model and should always check the restoration of the 0(3) 
symmetry. Unfortunately this was not done in [13, 14] and we shall see below that this symmetry 
is actually broken in the lattice models constructed by Catterall and Ghadab. 
The formulation in [13, 14] is based on a Q-exact deformation of the 0(3) sigma model in the 
unconstrained CP^ formulation. The authors used Wilson fermions which break supersymmetry 
softly and showed that the Q-based supersymmetry is restored in the continuum limit by studying 
the associated Ward identities. After performing a Hubbard-Stratonovich transformation to climate 
four-fermi terms the model contains two complex scalars u, a and Dirac fermions The path 
integral is 

•5b = 4 E {plid'^'^uUd'y^u), + pUAuUAu), + , (3.4) 

where the symmetric derivative in direction jj. is the arithmetic mean of the forward and backward 
derivates, 9^^"" = |(5+ + d~). The derivative operators without lower index denote the complex 
lattice derivatives d = di — 182 and B = di + 182- The lattice Laplacian A = d'^d~ originates 
from the Wilson mass term in the fermionic action 

Sf = ^^M[u,a]'^ (3.5) 
9^ 



with fermion matrix 



2 I iA - pu{Au) + h.c. ^^^-m _ 2pu(d'y'"u) + a 



^[^' ^] P ^Qsym ^ 2pu{d'y'^u) -a \A- pu{Au) + h.c. j ■ ^^'^^ 

Periodic boundary conditions ai^e assumed for all fields such that supersymmetry is not broken by 
the boundary conditions. The difference to a straightforward discretization is given by an improve- 
ment term 

AS = 4 E (^r^2 d'rui - d^ui dru2) q.d 

^ X 

which corresponds to the topological winding number and becomes irrelevant in the naive con- 
tinuum limit, similar to the improvement term in the M = 2 Wess-Zumino model [15]. For the 
simulation of this model on smaller lattices the HMC algorithm with an explicit calculation of the 
fermionic determinant is used. This has the advantage that no instabilities arise from introducing 
pseudo-fermions. Such instabilities may hide potential shortcomings of the lattice formulation. 
The improvement term is analyzed for a lattice size of = 8 x 8 at coupling 5"^ = 1.5. In our 
simulations a value of 5b ~ 2N is found (with statistical fluctuations of about 10%) in agreement 
with the simplest Wai^d identity [14], which is a consequence of the (nearly, up to the Wilson mass) 
lattice supersymmetry. In Fig. 6 the MC histories of and AS are plotted and they shed light on 
the influence of the improvement term on the dynamics. At a certain point in the simulation the 
value of ^3 freezes out and the improvement term starts growing largely negative. Just as for the 
Nicolai improved Wess-Zumino model [15] the lattice system is driven away from the continuum 
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Figure 6. MC history of (left panel) and AS (right panel) for a simulation of the lattice model (3.3) at 
coupling g^^ = 1.5 on an 8 x 8 lattice. 



physics, where AS" must vanish, into an unphysical phase . We conclude that the known problems 
of improved actions must be taken into account also in simulations of supersymmetric sigma mod- 
els - only configurations from the physical phase with nearly vanishing improvement term should 
be taken into account in measurements. Tunnel events to the unphysical phase are suppressed on 
large lattices and weak coupling g, i.e. in the continuum limit, similarly as for Wess-Zumino mod- 
els. Nevertheless, these observations suggest that similar problems may arise in simulations of 
all lattice models with exact supersymmetry that are constructed from a Q-exact action, e.g. two 
dimensional M = 2 Super Yang-Mills [28]. 

But why did this instability not show up in the results of [14]? The answer may be that in the 
rational HMC algorithm used in the simulations spectral bounds must be chosen to cover the spec- 
trum of M^M. Typically these are obtained by test runs with rather pessimistic bounds and small 
statistics, such that only the physical phase with /S.S w is present. But for the simulation that 
is shown in Fig. 6 the lowest eigenvalue of M'^ M decreases by a factor of 10^^ when entering 
the unphysical phase.^ For that reason the rational hybrid Monte-Carlo algorithm with spectral 
bounds that are not applicable to the whole simulation will be an inexact algorithm and will give 
an arbitrarily small acceptance rate for the unphysical configurations that dominate the path inte- 
gral. Furthermore sign det M is not positive definite and a deflated rational hybrid Monte-Carlo 
algorithm is necessary to get reliable expectation values. 

These results imply that for a reliable measurement of (^3) large g~'^ must be used in order to 
suppress the unphysical phase. Since the continuum hmit is reached for 5 — )■ 0, measurements will 
be affected by finite size corrections. The corrections of the observable n^, are assumed to vanish 
exponentially with growing volume, 

(7-^3)(iVs) = (n3)(oo) + Ae-^^% (3.8) 

such that a fit to this functional form for G {10, 11, 12, 13, 14, 16} and corresponding lattice 
volumes N = reveals the infinite volume value for g'"^ G {3.5, 4.0, 4.5, 5.0}, see Fig. 7 (left 
panel). Although (n^) decreases for fixed N^; and growing g'"^, the infinite volume values tend to 

'The normalized improvement term AS/N fluctuates around zero with a width of about 0.002 in the physical phase. 
^The largest eigenvalue of Al'' M is kept at the same order of magnitude in the unphysical phase. 
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grow for larger g ^, see Fig. 7 (right panel). Therefore the 0(3) symmetry will not be restored 
in the infinite volume continuum limit of the lattice model (3.3). All these results have a crucial 
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Figure 7. Left panel: (n^) for different lattice volumes N = at four couplings g~^. The infinite volume 
extrapolation according to Eq. (3.8) is indicated by the shaded areas. Right panel: Infinite volume values of 
(na) for different g~^. 

implication. Although the formulation (3.3) may restore the full M = (2, 2) supersymmetry on the 
lattice, the resulting continuum model is not the supersymmetric 0(3) model, because the global 
0(3) symmetry, that cannot be broken spontaneously in the continuum model [27], is not restored 
in the continuum limit. 

3.2 0(3) symmetric formulations 

The previous studies have stressed the importance of an 0(3) invariant formulation of the the- 
ory. The discretizations given in section 2.3 respect the global flavor symmetry and in order to 
check whether the coiTcsponding algorithms respect it as well we record the expectation value of 
n. For the coset construction with group valued field R the expectation value vanishes and the 
autocorrelation times are small. For the stereographically projected fields u we generate 0(3) 
symmetric configurations but observe large autocorrelation times for observables which depend on 
the field component coiTcsponding to the projection axis. This stems from the interplay between 
stereographic projection and molecular dynamics steps which introduce two pseudo-momenta and 
a pseudo-Hamiltonian for u to generate test configurations. Computing the pseudo-momenta for 
the constrained variables, we see that the momentum corresponding to the projection axis used in 
the HMC algorithm takes values roughly one half of the other momenta. In the language of uncon- 
strained variables, evolution of configurations is glued to a hyperplane of constant u^, which can 
be seen from figure 8 (left panel), where it^, v'^ and w'^ correspond to the three possible (canoni- 
cal) projection axes, while the actual projection axis of the HMC algorithm is fixed. To circumvent 
this problem, the projection axis in the molecular dynamics step is changed repeatedly between 
successive trial and acceptance steps. While it is in principle possible to choose random projection 
axes for every update, we rotate through the three canonical axes, saving only every third config- 
uration which corresponds to a fixed projection axis. A Metropolis acceptance step is yet needed 
after each change of the projection axis. The right panel of figure 8 illustrates how this improved 
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Figure 8. MC-histories for the lattice averages of u^, and using the three different projection axes 
in the usual HMC algorithm (left panel) and the improved algorithm (right panel). Note that fluctuations for 
which corresponds to the projection axis chosen in the HMC algorithm are severely suppressed in the left 
panel. 



update scheme restores the balance between the three projection axes. The configuration space is 
traversed quickly and the expectation values decrease with ongoing MC time. 

3.3 Fermionic masses 

The classical theory has no intrinsic mass scale. But there is a relation between mass gap and bare 
coupling by dimensional transmutation. In the MS scheme it is possible to compute the mass gap 
[8, 29] in relation to A^. Investigations of the M = 1 Wess-Zumino model have revealed that 
the fermionic mass is less affected by finite size effects than the bosonic mass. For that reason the 
fermionic mass is used to set the physical scale of the lattice, i.e. the lattice spacing and the physical 
box length. The 0(3) symmetric fermionic correlator for group-valued fields is constructed as 

(itpxtpy) = 2.-^ j VnVaVil) 5{nil))5{-n? - 1) {ixj^^ipy) e'^ 

= j VRVae"^^ j Vx iXx^'^ RlRy^^Xye-^' (39) 




with Sf given in (2.21) and 'trt^' indicates the trace over flavor and spinor indices. The coitc- 
sponding timeslice correlator is given by Cp{t) = N^^ J2xy {^'^{t,x)'^{o,y))- order to measure 
the mass in one of the ground states the configurations are projected, without loss of generality, 
onto the sector with H > for aE = J2x This can be achieved by flipping the sign 

of a for configurations with H < 0. Using these definitions the fermionic masses are obtained by 
a cosh fit to the correlator over the range t s]0, Ni[. A similar procedure is utilized in the case of 
stereographically projected fermions, where we have to use the projected fermion correlator instead 
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of (3.9), 



+ l6pl{Uxi\x){UxUy){XyUy)pl) , 

(itpx^lpy,!) ={lQpi{Ux\\x){\Uy)p^y) . 



(3.10) 



In this case non-diagonal elements contribute to the full correlator, which relate unconstrained 
fields of different flavor. A comparison of the fermionic masses for different lattice sizes using the 
SLAC derivative (see Figure 9) reveals that finite size effects on mpa are within statistical error 
bars if mpL > 5. 
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Figure 9. Fermionic masses from SLAC ensembles in units of the box length (left panel) and lattice spacing 
(right panel) for different lattice volumes and bare couplings computed using up to 2 • 10^ configurations. 



3.4 Chiral symmetry breaking 

Due to the orthogonality of Rx the chiral condensate in the coset formulation simplifies to 

(i^^V-x) = VRVae-^^ j VxiXxXxe-'^' = {trUCQ^^)) . (3.11) 

Using stereographic projection, however, we have to compute the projected condensate given in 
terms of unconstrained fields. This is obtained by using the trace of the correlator given in equation 
(3.10). If we rewrite the action (2.17) as S = SB[u,a] + X^PX, we can replace the quadratic 
fermion operator by 

{■■■>^^Ay,-) = {-C{P-')^y,y:)- (3.12) 

The continuum model is invariant under a discrete chiral symmetry i/; — )• ^7*1/^. It is sponta- 
neously broken in the infinite volume Mmit and the supersymmetric ground states correspond to the 
two ground states of this broken symmetry [30]. A discretization based on the SLAC derivative 
maintains chiral symmetry on the lattice by the cost of having a non-local derivative. For every 
finite lattice volume the expectation value (irptj^) will vanish and is hence not the appropriate mea- 
sure to trace a broken symmetry. One should instead analyze the histograms of the volume average 
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Figure 10. Left panel: Probability density of the volume averaged chiral condensate for lattice size 9x9 
at a bare coupling g^^ = 1 in the sign quenched and reweighted ensemble using the SLAC derivative. 
Right panel: Constraint effective potential (normalised to min(a2) U{aE) ~ 0) of the chiral condensate for 
different lattice volumes at g'^ ~ 0.64 computed using up to 3 • 10*^ configurations and the SLAC derivative. 

aS = ti^xfpx- Fig- 10 (left panel) clearly shows a double peak structure of the correspond- 

ing distribution p{aE), coinciding with the two ground states. The reweighting process reveals that 
a cancellation between positive and negative Pfaffians happens mostly for aH « 0. In the analysis 
of the constrained effective potential U{aE) = — ln(/j(aH)) /N for several lattice volumes at fixed 
coupling (see Fig. 10, right panel) no running of the two minima is visible, such that there will be 
a spontaneous chiral symmetry breaking in the infinite volume limit of the lattice model. 

3.5 Bosonic masses 

In order to test whether the chosen discretization corresponds to a supersymmetric theory, one 
should check the degeneracy of the bosonic and fermionic masses. The bosonic masses m^L are 
extracted from the 0(3) invariant correlator (2.31) via a cosh fit over the range t g]0, A''t[. The 
bosonic correlator is unaffected by a change in a, such that no projection on one of the two ground 
states is necessary. Having computed the boson masses, we can compare them with the fermionic 
ones for different couplings and lattice sizes. 

Simulations using SLAC fermions Calculations have been performed on lattice sizes N € 
{5^,7^,9^} over a coupling range g~'^ G [0.4, 1.2]. The direct comparison is shown in Fig. 11 
(left panel) and the results seem to be disappointing at first sight. The bosonic masses lie consid- 
erably below the fermionic partners and this deviation becomes even more pronounced for larger 
lattices. However, this does not necessarily imply that supersymmetry will be broken in the con- 
tinuum limit. Already for the simple lattice J\f=2 Wess-Zumino model with spontaneously broken 
Z2 symmetry and one exact supersymmetry the masses split in the strong coupling regime at finite 
physical box sizes [15]. E.g. for couplings for which the one-loop perturbation theory fails and for 
a box size mpL 10, a 20% splitting with a smaller bosonic mass is observed. From that point 
of view the supersymmetric 0(3) nonlinear sigma model could be similar to a strongly coupled 
J\f = 2 Wess-Zumino model. The finite size effects may be even larger and a mass splitting of 
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much more than 20% would not be surprising for mpL < 10. Only an analysis of the mass ratio 
mB/mp in the large volume limit can uncover a restoration of degenerate masses. This is exem- 
plarily shown for the results on the 5x5 lattice in Fig. 1 1 (right panel). Despite the fact that lattice 
artefacts are sizeable the basic mechanism becomes cleai\ In the limit of lai^ge volumes a relation 

JTlB , AM 
> = 1 



m^L = mpL — AM 



(3.13) 



with constant AM is found, such that the ratio tends to 1 and the masses will be degenerate in the 
infinite volume limit. ^ 
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Figure 11. Left panel: Direct comparison of bosonic and fermionic mass from the SLAC ensemble in units 
of the box size for three different lattice sizes. The dotted line denotes the case mp = mg. Right panel: 
Difference ttibL — mpL for varying box size nipL on a 5 x 5 lattice. The shaded area denotes a fit according 
to Eq. (3.13) for nipL > 6. 



But the accessible physical volumes at larger lattices still do not allow for a reliable extrapolation 
of the corresponding AM and no continuum limit can be taken at the moment. It is hence an 
open question if eq. (3.13) also holds true in the continuum limit and if supersymmetry will be re- 
stored. These questions should be resolved by computations on larger lattices. Such computations, 
however, become unfeasible because of the sign problem, which worsens with increasing lattice 
size in case of the SLAC derivative, see appendix D. Nevertheless, further information about the 
supersymmetric features of the implementation can be obtained by studying a corresponding Ward 
identity, cf. section 3.7. 

Simulations using Wilson fermions While the applicability of the SLAC derivative is confined 
to small lattice volumes due to the strong sign problem, Wilson fermions offer the possibility to 
explore larger volumes by utilizing efficiently preconditioned pseudofermion algorithms (compai^e 
to appendix E), however at the cost of larger lattice artifacts. Bosonic masses ai^e extracted as 
explained beforehand by calculating the stereographically projected correlators 

= 4:{PxUxUyPy), 

{nx,iny,i) = {px{l - ul) (1 - ul)py). (3.14) 



'For the 5 X 5 lattice a fit to Eq. (3.13) for mvL > 6 gives AM = 2.56(10). 
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For lattice volumes up to 32^ we obtain a dis- 
crepancy between bosonic and fermionic masses 
with lighter bosons, see Figure 12. This gap in- 
creases for larger lattices and there is no indication 
of a degeneracy in the continuum limit, so that su- 
persymmetry is seemingly not restored. This is 
not too surprising, though, since Wilson fermions 
may break supersymmetry in a way that is not dis- 
solved in the continuum limit. 
In section 2.5 we showed that it is not possible 
to construct a discretized action which simulta- 
neously respects at least one exact supersymme- 
try as well as the 0(3) symmetry. An improved 
action thus relies on adding fine tuning terms as 
a compensation for renormalized couplings that 
arise from symmetry breaking terms on the lattice. 

3.6 Fine tuning of the Wilson derivative 

From J\f=l Super Yang-Mills theories it is known that the correct continuum limit may be achieved 
by a fine tuning term that resembles an explicit fermionic mass such that the renormalized gluino 
mass is zero [31]. A similar procedure will be performed in the present case by deforming the 
fermionic derivative using a fine tuning mass m, 

= 7^(9^^'"),, + ^"'^ yA,, + 5"Pm5,y, (3.15) 

which enters the Hopping parameter k = (4 + 2m)~^. Additional degrees of freedom in the fine 
tuning procedure increase the numerical effort considerably, such that a careful choice of tuning 
parameters is necessary in order to keep the RHMC algorithm exact and efficient. In particular-, for 
different k, we have monitored the sign of the Pfaffian determinant as well as the spectrum of the 
Dirac operator, which is approximated by rational functions [32]. To guide our efforts, we measure 
the chiral condensate, bosonic and fermionic masses as well as the bosonic action for several values 
of the fine tuning parameter. 

The chiral condensate is extracted from the trace of the projected correlator, as described in (3.10), 
while the sign of the Pfaffian is taken into account by a reweighting procedure, 

(sgnPf 

V'V' = ^7 ^TT^. (3.16) 

(sgnPf)^ 

Here, denotes the sign-quenched ensemble. For the Wilson derivative we expect that one 
of the ground state energies is raised due to the explicit breaking of chiral symmetry, so that we 
get fluctuations around one of the minima in Figure 10 only. Switching on the fine tuning mass, 
the condensate is driven to larger values, showing a jump at some distinct point (see Figure 13, 
left panel). To gain a better understanding of this behaviour, we utilize histograms of the chiral 
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Figure 12. Comparison of bosonic and fermionic 
mass in units of the box size for three different lat- 
tice sizes using Wilson fermions. The dotted line 
denotes the case mp = "^b- 
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Figure 13. Left Panel: Expectation value of the chiral condensate aS for lattice volume 16^ and different 
couplings g~^. A jump corresponding to a first order transition is visible. Right Panel: Histograms of the 
chiral condensate aS for different values of the Hopping parameter (N ~ 24^, g^^ = 2). 



condensate to measure the distribution function p{aE), which is formally obtained by introducing 
a delta function into the partition sum, 

If 1 _ 

p{aE) = - VnVtP 6{aE - i^t/;) e"^ ^ 77 5^ -^(^^ " i(^V')i) • 

i=i 

It is henceforth possible to express the expectation value for the chiral condensate using this quan- 
tity as (ii/^tp) = X^i^i oSi/o(aSi)/ X^i^i /o(oSi)- Simulating the sign-quenched ensemble, we 
need to use the reweighted distribution function. 



1 ^ 



i=l 



sgnPft 
(sgnPf), 



(3.17) 



which may no longer be interpreted as a probability distribution, since sgn Pf can be negative. 
In contrast to the case of SLAC fermions, where 
no problem arises due to intact Z2-symmetry, the 
Wilson prescription leads to an additive renormal- 
ization, such that configurations with negative sign 
yield entries corresponding to a sign-flipped renor- 
malization constant. To avoid this behaviour, we 
will omit the sign information when considering 
histograms, seeing that the Wilson ensemble shows 
frequent sign fluctuations only in the direct vicin- 
ity of the critical Hopping pai^ameter and therefore 
allows to extrapolate from a region where omitting 
the sign is safe. Expectation values are however al- 
ways determined using the reweighting procedure 
and are hence not affected by this approximation. 
By an appropriate choice of the tuning parameter 
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Figure 14. Critical value of the Hopping parame- 
ter for several lattice sizes. 
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we obtain the spontaneously broken signature that is expected in the continuum limit (see Figure 
13, right panel), modified by the additive renormalization. The point of steep increase is identified 
as a first order phase transition in the chiral condensate, which is analogous to the case of Super- 
Yang-Mills theories. Using this signature, we have determined the critical value of the fine tuning 
parameter k for lattice sizes 8^, 16^ and 24^ and coupling g'"^ = 1 ... 2, see Figure 14. 
Regarding the masses of the elementary excitations, we see that the bosonic correlator is not af- 
fected by the fine tuning procedure and the bosonic mass takes a constant value within error bars. 
For the fermionic mass, we observe linear scaling behaviour for k < Kc and k > Kc, however 
in the vicinity of the scaling breaks down (Figure 15, left panel). Therefore, a reliable 

extrapolation of the fermionic mass based on small values of k is not possible, at least for the 
lattice volumes considered. We will hence utilize the bosonic mass to fix the physical box size 
uibL for our investigations involving fine tuning. In contrast to e.g. M=l Super Yang-Mills [31], 
simulations at the critical point are feasible since the theory inhibits a finite mass gap even in the 
continuum. In the vicinity of the phase transition, mixing between the two ground states occurs and 




tubL 



0.45 




Figure 15. Left Panel: Scaling behaviour of the bosonic and fermionic mass for N = 16^ and = 1.4 
with Kc — 0.382(1), marked by the black dot. The dashed line denotes the lattice cutoff of 1/16. Right 
Panel: Comparison of bosonic and fermionic masses in units of the box size for three different lattice sizes 



at K 



the itfjip correlator must be projected onto one of the ground states, analogously to the SLAC case 
(see chapter 3.3). However, since we may only determine the additive renormalization constant 
from the chiral condensate distribution function up to finite precision, fermion mass extraction is 
affected by a systematic error which comes from configurations with close to vanishing renormal- 
ized chiral condensate that are erroneously weighted as belonging to the false ground state. This 
error is negligible for large lattices and g^^ though, since these configurations are suppressed by 
the finite tunneling probability between both ground states. We observe a great improvement re- 
gai^ding the degeneracy of the masses (Figure 15, right panel), which is expected in the continuum 
Unlit from supersymmetry restoration. In particular, the fermionic masses no longer "run away" 
if the volume is increased, which hints at a proper cancellation of the divergent operator causing 
these problems (compare to Figure 12). Nevertheless, a true proof of this conjecture may only be 
provided by a study of all divergent operators based on lattice perturbation theory, which is not 
pursued here. For finite box sizes of mj?L > 2 a thermal mass-splitting similar- to the SLAC case 
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seems to emerge, however showing a shghtly smaller gap. To explore this region, further large 
volume simulations would be needed in order to suppress lattice artifacts. 



3.7 Path integral based Ward identity 

Similar- to the Ward identity that is given in [14] an equivalent relation can be derived for the present 
lattice model. Starting from the path integral in the continuum, Q exactness is given by the twisted 
supercharge [14] in the continuum, such that S = ^QA with = 0. This implies the continuum 
Ward identity 

|^ = (-iQA> = 0, (3.18) 

since action and measure are invariant under the symmetry transformation generated by Q. After 
putting the theory on the lattice, integrating out the constrained auxiliary field /, and introducing 
the unconstrained field a to get rid of the four-fermion interaction (in that order) one is left with 
the path integral, 

Z = g^ j Di? Do- DA e-^[^''^'^l (3. 19) 

with as number of lattice sites and the action given in eq. (2.20). Here, the coupling dependent 
part is important, so that only constant numerical factors may be dropped. The derivative of the 
Schwinger functional is 

where Qp denotes the Dirac operator in (2.22), and dim is its dimension in terms of flavor and 
spinor components as well as lattices sites. In our case the Ward identity (3.18) reads 

(^b) = (3.21) 

with 5b defined in Eq. (2.21). The same Ward identity can be derived for the formulation based on 
the stereographic projection. 

In order to see a possible restoration of supersymmetry in the continuum limit, the bosonic action 
has been calculated with the SLAC derivative for three lattice sizes in a coupling range where finite 
size effects should be negligible, i.e. for m^L > 5. The results that are shown in Fig. 16 (left panel) 
reveal that for the smallest (5 x 5) lattice Eq. (3.21) is violated as much as 10% for mpL ^ 5 and 
up to 20% at mpL 10. Therefore the Ward identity violation grows for coarser lattice spacings. 
However, in the continuum limit at a fixed physical volume m^L the Ward identity tends to be re- 
stored, cf. Fig 16 (left panel). Additionally one can explore the small volume regime of this theory 
by sending ^ — )• at fixed lattice volume to reach the continuum limit. This has been performed 
on a 7 X 7 lattice for a large range of couplings e [0.4, 100], see Fig. 16 (right panel). Here 
the Ward identity is explicitly restored in the limit of large g~^, although this result has to be taken 
with care since the physical box size becomes unreliably small. 

The examined lattice sizes are rather small, but these observations at least imply that a supersym- 
metric continuum limit can be reached and that the non-degeneracy of bosonic and fermionic mass 



The factor g in front of the path integral stems from the Gaussian integrals that need to be carried out for / and be 
introduced for a. 
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Figure 16. Normalized action 2 (5b) /3A^, which will take the value 1 if supersymmetry is restored in the 
continuum limit as required by the Ward identity (3.21). Left panel: Measurements for different lattice sizes 
with physical volumes 5 < mpL < 1 1 in the SLAC ensemble. Right panel: Results for fixed lattice volume 
TV = 7^ (SLAC) and iV = 8^ (Wilson, k = 0.25 fixed) and couphngs g-'^ € [0.4, 100] that reach the regime 
of small physical volumes (small g). 



is a finite size effect. However, in order to obtain a definite answer simulations on larger lattices 
would be necessary. For instance, one ought to verify that (5b) /N does not undershoot and drop 
below 1.5 in the continuum limit at fixed mpL. Using the Wilson derivative, we measured the 
bosonic action in the finetuned ensemble at k = Kc and observe similar behavior, although for the 
largest lattice considered we see a discrepancy of 7% for fine lattice spacing (g'"^ = 2) and up to 
14% for a coarser lattice spacing (g""^ = 1.4), see figure 17. However, unlike for SLAC fermions, 
the Ward identity approaches its continuum value from above in the limit of to infinity and 
therefore is unlikely to undershoot (see Fig. 16, right panel). Hence a comparison of both deriva- 
tives is only valid in the regime of large g'"^ where monotonic behaviour is encountered, which is 
however unavailable for the SLAC derivative due to either fine-size effects for small lattices or the 
strong sign problem for large lattices. Overall, the slope of the Ward identity clearly points to a 
restoration of supersymmetry in the continuum limit. 
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Figure 17. Expectation value of the bosonic action at k = Kc for different g ^ (left panel) and different box 



sizes (right panel). 
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1 is expected in the supersymmetric continuum limit. 
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4 Conclusion 



The purpose of this article was to obtain a lattice formulation of the supersymmetric nonlinear 0(3) 
model which maintains the symmetries of the theory, at least in the continuum limit. The tai^get 
manifold of the model is Kahler and hence it possesses an = 2 supersymmetry, which could 
provide for a nilpotent supercharge and a discretization prescription based on this which maintains 
an exact supersymmetry on the lattice. We first derived the supersymmetry transformation in terms 
of constrained field variables and used these to analyze the applicability of such an approach in the 
nonlinear 0(3) model. We realized that there is in fact no way to formulate a discretization of the 
theory which maintains simultanously both the 0(3) symmetry and a part of supersymmetry. We 
explicitly verified that the Q-exact supersymmetric formulation presented in [14] breaks the 0(3) 
symmetry such that it is not restored even in the continuum limit. 

In contrast, we started from a manifestly 0(3) symmetric discretization and investigated whether 
the supersymmetry is restored in the continuum limit. The spherical geometry of the target space 
is treated by two separate approaches, group valued variables on the one hand and a stereographic 
projection on the other hand. The numerical simulations aie performed with the SLAC derivative 
and the Wilson derivative respectively. The former does not break the chiral Z2 symmetry of the 
classical action explicitly, allowing for the evaluation of histograms to verify the symmetric ground 
state structure which is spontaneously broken. The applicability of the SLAC derivative in theories 
with curved target space was illustrated by a calculation of the step scaling function in the quenched 
model. A disadvantage of the SLAC derivative, however, is the strong sign problem that becomes 
relevant already at comparably small lattices. In a sense the SLAC derivative is already too close to 
the continuum limit since it conectly reflects the intrinsic sign problem even on moderately sized 
lattices. In order to test the supersymmetric properties of our lattice models we analyzed the masses 
of fermions and bosons as well as a Ward identity based on the bosonic action. For both derivatives 
the results indicate that the Ward identity is fulfilled in the continuum limit at finite (large) physical 
volume. Concerning the expected degeneracy of the masses, no final statement can be done based 
on the small lattice sizes that are accessible by the SLAC derivative. The Wilson derivative en- 
ables us to investigate larger lattices, but breaks chiral symmetry explicitly at finite lattice spacing, 
which leads to the renormalization of relevant operators in such a way that one is carried away 
from the supersymmetric continuum limit. Motivated by Super Yang-Mills theory, we fine tuned 
the fermionic mass in order to remove the explicit breaking of the chiral symmetry. Furthermore, 
using the thus tuned ensemble, no "run away" of fermionic masses is visible, but rather a degen- 
eracy of boson and fermion masses. This indicates that a single parameter is sufficient to provide 
for a supersymmetric continuum limit and for a cancellation of the encountered divergences. This, 
however, remains a conjecture until explicitly checked by means of lattice perturbation theory. 
We have therefore presented a lattice discretization which incorporates the 0(3) symmetry exactly 
at finite lattice spacing and furthermore shows restored supersymmetry as well as the spontaneously 
broken chiral symmetry in the continuum limit. The price to pay is a single additional fine tuning 
parameter It remains an open question whether further possibilities for the lattice derivative like 
the Neuberger operator provide a discretization which is free from explicit chiral symmetry break- 
ing and perhaps free from the need for fine tuning, provided that only a mild sign problem is 
encountered at the same time. 
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A Conventions and Fierz identities 

We choose the Majorana representation 

70 = 0-3, 71 = -f^i; 7* = 17071 = cr2, C = -i(T2, (A.l) 

and the conjugate spinor is defined as x = 'X^C. In the main body of the paper we employ the 
Fierz relation 

V'X = -kx^- i(X7^^)7m - \{X1*^)1*- (A.2) 
Due to the symmetry properties 

xi^ = iPx, xi^i^ = -i'l^x, xi*'4' = -V'7*x (A.3) 

the two last terms in (A.2) vanish for x = such that 

ipip = -^'ijjipl . (A.4) 

B Invariance of the action under the second supersymmetry 

We will prove the invariance of the on-shell action 

S[n,ip] = j (fx {d^n • + i-0^^ + i (tp^)^) (B.l) 

under the second supersymmetry transformations (2.12). The variation of the Lagrangian is^' 

6C = 2\d^n ■ a^(n x eV) - 2i'4)${n x $n£) + 2ip0{eip x tp) - (tptp) 'ip{n x 0ne) . (B.2) 

The term oc -0^ vanishes, since -0 is a Grassmannian field with only four independent degrees of 
freedom. We will see that the first and second term cancel each other as well as the third and fourth 
one. Starting with the first two terms, they can be written as 

2i9^n(n x ed'^ip) - 21-0(71 x d^ne) - 2['ipY'Y£{df,n x (9^n) . (B.3) 

The last term vanishes since d^n x d^n is parallel to n and hence perpendicular to ip. Integrating 
the second term by parts one sees that the first and second term cancel owing to the cyclicity of the 
triple product. 

"up to a negligible boundary term i/)7^(£i/) x V")) 
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The cancellation of the third and fourth term in (B.2) is a bit more involved. First, we partially inte- 
grate the third term and obtain — 25^i/j7'^(ei/; x -i/?). Since n • •0" = for both spinor components 
a, we conclude that etf^ x ip is parallel to n such that 

- 25^'07^(e'0 xip) = -2(9^'07''n) n{eip x tp) . (B.4) 

The condition xpn = implies dfj^ipn = —xpd^n. (B.4) can hence be written as 2{ip0n) n{eip x 
ip). To proceed further, we utilize V^irii = —^^2^2 — 'ipz'n^ and the Fierz relation Tpi'y'^ipi = 0: 

2{ijji0ni + ■02^f^2 + ip30ri3)[nie'4>2 ■ V'3 - niSip3 ■ ^2 + cyclic terms] (B.5) 
= 2^/527'' V's (9^ "-2 • ni - d^ni • 712) eV2 + 2^37^^2(5^?^l • - 5^?^3 • ?^i)e^3 + cyclic terms . 

Finally, we employ the Fierz relation (07^/3) ea = ^aa{j5^^e), which holds for Majorana 
spinors, and obtain 

{4'2'4'2)i''il^e{nid^n2 - n2d^ni) + (^3V'3)V'27^'^("'3<9^"-i - nid^n^) + cyclic terms . 

Finally, using {aa)a = we obtain for the third term in (B.2) the simple expression 

{xpip) ip{n X $ne) . (B.6) 

As a result, the third and fourth term in (B.2) cancel each other. This proves that the action is 
invariant under the second supersymmetry transformation (2.12). 
The invariance of the constraints can be shown easily: 

5{n^) = 2m ■ {n x eip) = 
6{n • tp) = i{n x eip) ■ ip — n ■ {n x $ne) — in • {e%p x tp) = . 



C Transformation of the discretized measure 

The transformation from the constrained fields (n, ■?/>) to the unconstrained fields {u, A) has a non- 
trivial Jacobian. Since the transformation only relates values of the fields on a fixed lattice site it is 
sufficient to calculate the Jacobian for a given site. Denoting values of the fields on this site by 

we are lead to consider 

6{n^ - l)5{n ■ ip^)5{n ■ ip^) 



2\ni 



ni 

Consequently, the measure on a given site transforms as 

dndip^dtp^ 5(n2 - l)6{n ■ ip^)6{n • -0^) = i J(ii) dit dAMA^ (C.l) 
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with Jacobian 



J(-it) = Jl-n\{u) |sdet{(n_L,'0_L) ^ (w,A)}| . (C.2) 



In the super-stereographic projection (2.7) n does not depend on A°, and -0° does not depend on 
A'^ for [i ^ a. The superdetemiinant is hence given by 

.det{(n,.^,) (..A)) = ;gf,"-f (t.^,,.) ■ (C3) 
In an O(A^) model, all three determinants are equal to 



{2pr-^—— with p 



1 + 1*2 1 + ^2 

Expressing the square root in (C.2) in terms of the new fields, 



1-1*2 

we end up with the Jacobian 

JN=(^;^°^(i + «')'^"'- (C-4) 

The functional integral measure for the supersymmetric 0(3) model with uncontrained fields is 
thus 

Wdu^dXldXl {l + ulf . (C.5) 

X 

Note that we proved on the way that the Jacobian for the purely bosonic O(A^) model is 

JbHoc/,^-i = ^^^^1^. (C.6) 



D Sign of the fermion determinant 

In order to check whether the sign-quenched approximation is applicable in case of the SLAC 
derivative, simulations on lattice sizes ranging from 5x5tollxll have been performed over 
a coupling range g~'^ G [0.4, 1.2]. The results that are based on 10^ configurations per data point 
(see Fig.18, left panel) indicate that the average sign of the Pfaffian is smaller for smaller g""^, 
which is equivalent to coarser lattices. The problem is that the sign problem worsens for larger 
lattice volumes at fixed coupling. In these cases the probability based, i.e. sign-quenched, Monte- 
Carlo sampling will not correspond to the relevant configurations in an unquenched ensemble and 
statistical errors on reweighted measurements will become rather large. Nevertheless, with standard 
Monte-Carlo techniques simulations are only possible without taking the sign into account, such 
that a reweighting becomes unavoidable. The discretization based on the SLAC derivative becomes 
thus unfeasible for larger lattices, which we require in order to study the continuum Umit and a 
possible restoration of supersymmetry there. 

As a consequence, we have to rely for this purpose on the Wilson derivative only. The average 
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Figure 18. Left Panel: Average sign of the SLAC Pfaffian for different couplings g^^ on lattices sizes 
ranging from 5 x 5 to 11 x 11. Right Panel: Sign of the Wilson Pfaffian on a iV = 16^ lattice for different 
couplings g^^ and normaUzed finetuning parameter k/kc- 



Pfaffian sign for tiie Wilson derivative is depicted in tlie right panel of Fig. 18. We see that 
for small values of the finetuning parameter k, which correspond to little finetuning, changes in 
the Pfaffian sign are suppressed and it is possible to evaluate expectation values directly without 
the need of reweighting. In the vicinity of however, this behaviour changes from a mild 

correction for large coupling g^^ = 2 up to a significant correction for smaller coupling = 1.2. 
In the finetuned ensemble at k = Kc, we regain the sign problem known from the SLAC derivative. 
Still, applying the Wilson formulation allows for an efficient even-odd preconditioning scheme, 
which provides a sufficient speed-up in order to account for the larger configuration numbers that 
are needed here. 

E Algorithmic aspects 

We argued in section 3.2 that it is of considerable importance to tune the many technical parameters 
of advanced algorithms like Rational Hybrid Monte Carlo (RHMC) in the right way in order to ob- 
tain an appropriate algorithm that parses configuration space in the correct way. With the sign prob- 
lem and lai^ge condition numbers lurking, a Hybrid Monte Carlo algorithm with exact evaluation 
of the fermion determinant and the inverse fermion matrix using efficient LAPACK routines was 
used at small lattice volumes up to 16^ in order to provide a solid ground for the implementation 
of more sophisticated pseudofermion algorithms, e.g. RHMC. These algorithms become necessary 
at large volumes due to the poor scaling behaviour of the LU decomposition used to solve for the 
inverse fermion matrix. However, no efficient preconditiong schemes exist for the nonlocal SLAC 
derivative. Therefore we concentrated our efforts on the Wilson derivative which allows for even- 
odd preconditioning and, moreover, shows less severe sign problems at finite volume. In Table 2 
average condition numbers k obtained from the exact matrix norm k = \\Q\\ • ||Q~^|| are shown 
as well as typical iteration numbers of a conjugate gradient solver for three different choices of the 
fermion matrix: 1. the original fermion matrix, 2. the reduced fermion matrix where a factor p, 
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original 


reduced 


even-odd 


condition number 


1.6(6) • 10^9 


1.3(6) • 10^ 


1.4(7) • 10^ 


eg solver steps 


454(10) 


152(1) 


48(1) 



Table 2. Average condition numbers and CG solver steps for three different choices of fermion matrix 
formulation. 



whose determinant can be evaluated analytically, is separated on both sides, 



Q' 



xy,ij 



4((5jfc - 2ux,iUx,kPx)M"l^{Skj - 2uy^kUy,jPy) 



4.(3aJ,y5ij5''^ 



(E.l) 



and 3. the preconditioned reduced fermion matrix using the well-known even-odd preconditioning 
scheme [33]. The critical step of the iterative eg solver is the application of the inverse of Q^Q, 
which is used in the pseudofermion algorithm, to some random (pseudofermion) vector, Y = 
(Q^Q)-iX. We see a significant improvement in the number of solver steps and in the condition 
numbers for the reduced fermion matrix. This can directly be verified by looking at the eigenvalue 
spectrum depicted in figure 19. Both matrices, the original and the reduced one, are real and 
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Figure 19. Eigenvalue frequency of the original and reduced fermion matrix (left panel, normalized such 
that both the largest eigenvalue and the integrated surface equal one) and for the even-odd preconditioned 
matrix (right panel, no normalization, logarithmic scale) from a sample of 2000 configurations on a 12^ 
lattice. 

antisymmetric, so that all eigenvalues are purely imaginary and come in complex conjugate pairs. 
Whereas the original fermion matrix exhibits a large number of eigenvalues very close to zero^^, 
this is not the case for the reduced one and the condition numbers hence decrease drastically. A 
further improvement is achieved by even-odd preconditioning. The preconditioned matrix is no 
longer antisymmetric and we see that the majority of eigenvalues lie close to 1. This is what we 
expected, since the goal of the incomplete LU preconditioning [34, 35], which can be regai^ded as 
a generalization of the even-odd scheme, is to rewrite the fermion matrix in the form 1 — L — U, 
where 1 is the identity matrix and L (U) is a lower (upper) triangular matrix. 



^This is expected according to the Banks-Casher-relation. 
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